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Abstract 

We consider logarithmic extensions of the correlation and response functions of 
scalar operators for the systems with aging as well as Schrodinger symmetry. Aging is 
known to be the simplest nonequilibrium phenomena, and its physical significances can 
be understood by the two-time correlation and response functions. Their logarithmic 
part is completely fixed by the bulk geometry in terms of the conformal weight of the 
dual operator and the dual particle number. 

Motivated by recent experimental realizations of Kardar-Parisi-Zhang universality 
class in growth phenomena and its subsequent theoretical extension to aging, we in- 
vestigate our two-time correlation functions out of equilibrium, which show several 
qualitatively different behaviors depending on the parameters in our theory. They ex- 
hibit either growing or aging, i.e. power-law decaying, behaviors for the entire range 
of our scaling time. Surprisingly, for some parameter ranges, they exhibit growing at 
early times as well as aging at later times. 
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1 Introduction 

Most of the physical systems are out of equilibrium PQ[2]- Thus it is very interesting to 
understand various physical properties of those systems. In general, it has been known to be 
difficult to understand the nonequilibrium phenomena despite active studies. Still, progresses 
have been made in the context of the simple nonequilibrium critical phenomena, where one 
can use their underlying local scaling invariance, universality classes and critical exponents to 
capture universal properties. One simple example is Kardar-Parisi-Zhang (KPZ) universality 
class [3]jl]. This universality class is known to describe many interesting systems in nature, 
such as paper wetting, epidemics and fluid flow of porous media and so on [1] [5]. In the last 
few years, this class is experimentally realized in the context of surface growth in turbulent 
liquid crystals [6J[7J, which provides the experimental evidence of the dynamical scaling of 
the KPZ theory in 1+1 dimensions. 

Inspired by these developments, the authors of [8] investigate the possibility of aging 
phenomena, which can arise in the same KPZ universality class. It is shown that the systems 
undergo simple aging in the correlation and response functions, whose dynamical scaling 
is characterized by some aging exponents. Moreover, it is claimed that the form of the 
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autoresponse scaling function is well described by the logarithmic extension of the scaling 
function with local scale invariance [9]. This provides us the main physical motivation 
to pursue the logarithmic extension of the aging correlation functions in the context of 
holography, following [10] [11]. In the rest of the introduction, we briefly explain aging and 
logarithmic conformal field theories (LCFT). 

Aging is one of the simplest nonequilibrium phenomena of many body systems, which is 
realized when the system is rapidly brought out of equilibrium. The basic properties of aging 
can be characterized by two-time correlation functions with two-time scales, "waiting time," 
when the system is perturbed after brought out of the equilibrium, and "response time," 
when the perturbation is measured. The correlation functions show the typical behaviors: 
older systems with a bigger waiting time relax in a slower manner than younger systems 
with smaller waiting time, see e.g. [2J[5J[12J for reviews. 

In the context of Anti-de Sitter space / Conformal field theory correspondence ( AdS / CFT) 
[13] . the geometric realization of aging is put forward in [10] by generalizing Schrodinger 
background [H][l5] with explicit time dependent termsQ The authors of [10] observe that 
aging algebra can be realized as a subset of Schrodinger algebra using a singular time depen- 
dent coordinate transformation. Thus it is established that the time translation symmetry 
is broken globally, and the aging symmetry is realized as conformal Schrodinger symmetry 
modulo time translation symmetry. Yet, it is claimed that the time itself necessarily has to 
be complex to show a relaxation process, which is one of the central properties of the aging 
[TU] . This complexification of the time is certainly not preferred. 

This drawback is overcome in [11J, by introducing some decay modes of the bulk scalar 
field along the 'internal' spectator direction £, which is not explicitly visible from the dual 
field theory with Schodinger and aging symmetry. The resulting two-time correlation func- 
tions show a dissipative behavior and exhibit the characteristic features of the aging system: 
power law decay, broken time translation invariance and dynamical scaling between the time 
and spatial directions [llj . see also e.g. [2] [5] [12] from the field theory point of view. These 
properties are shown in the geometric realizations of the aging background in the context 
of [H][15] as well as of the aging in light-cone in the context of [T7][TH]. The corresponding 
finite temperature generalizations for these two different backgrounds with asymptotic aging 
invariance are also studied in [TTj . 

LCFT is a CFT, in which some correlation functions have logarithmic singularities. For 
the review on LCFT, see, for example, [19j[20j. Recently, there have been interesting de- 
velopments in LCFT as a dual CFT of higher derivative gravity model. It started in three 
dimensional gravity models [21]- and then in their higher dimensional analogues, namely, 

1 See also [TB] for some general time dependent deformations of Schrodinger backgrounds. 
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critical gravity [26]-[30j. Generically, these theories include higher derivative terms as well as 
Einstein-Hilbert and cosmological constant terms and admit the AdS spacetime as a vacuum 
solution. In four and higher dimensional spacetime, the higher derivative terms are given 
by some combination of R 2 and R^R^ terms. When the coupling of the higher derivative 
terms is tuned, the linearized perturbations of the metric around the AdS vacuum solution 
include the massless and logarithmic modes [26] |^] This is called critical gravity and its dual 
CFT becomes LCFT which includes stress-energy tensor and its logarithmic counterpart. 

The holographic descriptions of LCFT have been attempted since early days of the 
AdS/CFT correspondence using a pair of scalar fields [34J [35J [36] . Recently, the correla- 
tion functions of a pair of scalar operators in non-relativistic LCFT have been considered in 
the context of AdS/CFT correspondence using the Lifshitz background [37!]. In this paper, 
we study the correlation functions of a pair of scalar operators, which would play the role of 
order parameters, in the non-relativistic LCFT with Schrodinger and aging invariance. 

We introduce a new type of action with Schrodinger invariance and generalize the correla- 
tion function to obtain the logarithmic corrections in detail in section ^2} This is generalized 
for the system with aging invariance in section £|3} In section £|4j we investigate various 
physical properties of the logarithmic aging correlation functions and also seek possibili- 
ties to connect to the physical systems with surface growth [6] [7] and their extension with 
logarithmic corrections [8]. We conclude in section $5] 



2 Schrodinger LCFT 



Let us consider some strongly coupled (2 + 1) dimensional dual field theory with Schrodinger 
symmetry, living at the boundary of the Schrodinger background 

L 2 

dsl = — ( dy\ + dy\ — 2dx + dx~ dx +2 + du 2 ) , (1) 

u 1 V u J 

where 7 indicates the fact that we consider two different cases on the same footing, the 
Schrodinger background with 7 7^ [14J[JS] and the AdS in light-cone 7 = [T7][TB] at 
zero temperature. It turns out that these two geometric realizations give the same physical 



2 In the theory without these fine-tuned coupling of the higher derivative terms, the fluctuations around 
AdS spacetime include the massive and massless graviton modes, one of which is ghost-like. Since AdS 
spacetime has timelike boundary, one may impose a boundary condition and truncate the ghost-like massive 
mode if it falls off more slowly than the massless one[3T][3U]. It was shown in [35] [33] that the theory with 
the truncation becomes effectively identical to the usual Einstein gravity at the classical but full nonlinear 
level and thus the correlation functions of the dual CFT can be readily read off from those in the CFT dual 
to Einstein gravity. If we choose special values on the coupling of higher derivative terms, the massive and 
massless modes become degenerate and turn into the massless and the logarithmic modes [25] . 
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properties for the correlation functions at zero temperature for the scalar operators of given 
conformal dimensions. Their finite temperature generalizations have been realized in [38J- 
[12] for 7 7^ and in [43j[44] for 7 = 0. It is also observed that there exist clear distinctions 
in metric fluctuations between these two backgrounds 



Motivated by the recent interests on LCFT from the holographic point of view [37], we 
consider two scalar fields and in this background, whose action has the following form 

S = K j d A x jf du^g [g MN «9 A /0 d N( p + m 2 + ^0 2 ^) , (2) 

where ub represents a cutoff near the boundary, K = — 7t 3 L 5 /4k 2 , M,N = +,—,u,yi,y2 
and //, v — +,yi,y2- Upon integrating out the field 0, one can obtain an equivalent action 
with second order time derivatives and fourth order spatial derivatives. In the context of 
AdS/CFT, the logarithmic extension on Lifshitz backgrounds can be found in [37], where 
one can find some general comments for this action. 

Let us comment the general solution space of the theory described by the action Q. 
In [37], a bigger solution space for the Lifshitz theory with two spin fields and with 
different parameters m and m is suggested. When specialized with m — m, the theory 
describes the Lifshitz LCFT. It is straight forward to follow this features and generalize the 
solution space of the theory by introducing two different parameters m and m for the fields 
and 0. 

Contrasted to Lifshitz theories, Schrodinger theories have additional features due to the 
isometry associated with the coordinate x~, called "dual particle mass." This can be shown 
explicitly as <9_0 = —A4<j) and <9_0 = —Ai<f) considering the fields and with two different 
values of Ai and Ai, respectively^] Thus the solution space of our fields and can be 
described by the parameters (m, Ai) and (m, Ai). In this paper, we consider the case m = fh 
and Ai = At, which gives the Schrodinger LCFT. 



The linearized field equations for and of the action rt2h become 



2.M<9 + = -^D0-V 2 0+^0, (3) 
2Md + 4> = -\vj>- V 2 0, (4) 



v? 



where 

V = u 2 



3 , jM 2 -m 2 L 2 

2 



(5) 



u u z 

3 As mentioned in the introduction, we are concerned with the aging properties which shows dissipative 
properties in the correlation functions. Thus we consider only the decaying modes of the scalar fields cj> and 
(j). For the detailed discussion for the differences between the decaying modes and oscillating modes, see 
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which is a differential operator for the Schrodinger case. 

Following [H][37], we construct the bulk to boundary Green's function Gij(u,u>, k) as 



d 2 k did 



-ik'„x» 



G r n(w,A; At )J(A; /Lt ) + G 12 (u,k^J(k^) 
G 21 (u,k')J(k') + G 22 (u,k')J(k') 



(6) 



(2tt) 2 2tt 
d 2 k' duj' 
(2vr) 2 2^ 

where we use x^ = (x + ,y), k^ = (-co, k), which give k^ = k-y — ux + . We choose G 2 \ = 0, 
which is in accord with the structure of the equations of motion given in Q. The Green's 
functions satisfy 

(7) 



VG n = , 



T>G\ 2 — G 22 , 



VG 



22 







where T> is defined in (5) and q = \/k 2 + 2Aiiu. The Green's functions Gu and G 22 are 
given by 

Gn(u, k^) = en u 2 K v (qu) , G 22 (u, k^) = c 22 u 2 K v (qu) , (8) 



where v = y/4 + m 2 L 2 — 7.M 2 . The normalization constants Cn = C22 = u i ^(qug) can ^ e 
determined by requiring that Gh(ub, k^) = G 22 (ub, k^) = 1 [35] . 

There exists another Green's function G12 due to the action we choose ([2]), which satisfies 

VG 12 = G 22 . (9) 

To evaluate Gi 2 , we use the same methods used in [35]. Using 

d 



dv 



2v 



and the fact that VG 22 = 0, we get 



V 



1 d 
2v dv 



G 



22 



G 



22 



Thus 



G -——G 1 d ( u2R ^ \ 
12 2v dv 22 2v dv \u 2 B K u {quB) J 



(10) 



111 



12) 



We will not write down the explicit expression, as it is lengthy and not so illuminating. 

After plugging the bulk equation of motion into the action (|2]), the boundary action 
becomes of the form 



S 



L 5 ~u 2 
u = -K / d 3 x— <pj^d u (, 



K / dx + 9(x-' 



duo duj' 
2tt 2tT 



-i(u)—uj')x + I j 



d 2 k d 2 k! 
(2vr) 2 (2tt) 2 



x J(k' fl )[J : 'i(uB,k' fit ,k tl )J(k fl ) + J r 2 (u B ,k' fl ,k II )J(k fj 



(13) 



II D 
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For the Schrodinger case, the system has local time translation invariance, thus the time 
integral is trivially evaluated to give delta function. The J-"'s are given by 

Tx = ^G 22 (^) d u G u (^) , (14) 
^2 = ^G 22 (k' IM ) d u G u (y . (15) 

T\ can be expressed, using 

d u [u 2 K v {qu)\ = u[{2 - u)K u {qu) - quK v _ x {qu)\ , (16) 

as 



L 3 [ K v _ x (qu B )~\ 2Y(l-u)L 3 rqun\ 2v 



(2 - v) - qu B : 



I» u% V 2 ) 



m • (17) 



K v (qu B ) 

Note that we keep only the first non-trivial contribution in the small u B expansion and 
use the normalization for the wave function to be normalized to unity at the boundary 
|45j . This result is already noticed in several different places including [15] [IT] [11] . The 
momentum space correlation function is given by 

(<j>(u/, k')<p(u, k)) = -2(2n)- 3 5(k' - k) 5(u/ - u) Ti(u B: u, k) , (18) 

where we use (p(x + ,y) as the dual operator of the bulk field (p(u,x + ,y). 

The coordinate space correlation function can be computed using the inverse Fourier 
transformation 

Y{l-v) L 3 M 1+ » 6(xl)9(xl - xf) ( M{y 2 -y x f \ 

' eX P I nf + +\ I ■ V 1 ^/ 



T(u)T(-u)-k2"- 1 u b - 2u (xZ-xi) 2 +» 2(x+-x+ 

This is the result we are looking for. Note that the final result has the same form as that of 
the Schrodinger case which is given in [TT] . 



Now let us turn to evaluate the T 2 using (12) 

(0(4, &)<Kxi, m)) = -29(xi) [ ^e^Me^)^,,,^ 
1 L 3 M l+v 6(x+)6(x+ - xf ) ( M(y 2 -y x ) 2 



T(u)n2^u B - 2u (4-x+) 2 +- ' GXP V 2{xt-xi 



x l--V4 + z/ln[ , , " , J . (20) 
V T{y) [ 2(x+-x+)V 1 ' 

This is our main result in this section. The correlation function of the dual field theory 
acquires the logarithmic time dependence, which is very interesting. We give comments on 
it after we get more general results in the aging case below. 
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3 Aging LCFT 



In the previous section, we evaluate the two point correlation function of scalar order param- 
eters of the dual field theory with Schrodinger symmetry. Similarly, we would like to compute 
the correlation functions with aging symmetry, which is intrinsically out of equilibrium and 
time dependent. 

Let us consider the following aging background 

ds 2 u = ^ (dyf + dy\ - 2dx + dx~ - (-^ + dx +2 + —dudx + + du 2 ] , (21) 

whose construction is explained in [I0][TT]. The parameter a signifies the fact that the 
system has an explicit time dependence, while the parameter 7 represents the difference 
between the aging in Schrodinger background with 7 = 1 and the aging in light-cone with 
7 = 0. The system is one of the simplest time dependent background, with which we can 
explicitly explore important physical properties related to out of equilibrium. The reader 
can find more detailed information in [TU] 

Consider two scalar fields and whose action is given by The linearized field 
equations for <fi and (f> become 



2M [d + + 
2M (d+ + 



aM 
2x~+ 
aM 
2x~+ 



1 - 9 1 ~ 

— -V a <j) - V 2 + 



u- 



--V a <j)-V 2 



where we treat x coordinate special by assigning d-4> = —M<t> and <9_ 



u 



dl- 



2Ma + Z n . 4Ma + a 2 M 2 + 7M 2 -m 2 L 21 
u u 2 



-d u + 



(22) 
(23) 

-M(f> and use 
(24) 



which is a modified differential operator for the aging case. 

The bulk to boundary Green's function GfAu, u>, k) may be defined similar to the Schrodinger 
case given in 

0(w,£ M ) = 



d 2 k du; 



[u, x h 



(2tt) 2 2tt 
d 2 k' du' 



— ik' 
e * C 



G^(u,k')J(k') + GUn,k')J(k[ 



(2vr) 2 2tt 

Once again, we choose Gfi — and the Green's functions satisfy 

, T) n G°tn — Goo 1 T^nGoo = 



VaG^ 



(25) 
(26) 

(27) 
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-aM 



where V a is defined in (24). We find that Gfj = ^^-J G^, where G^ satisfies equations 
g - Q of the Schrddm erer case. 

After plugging the bulk equation of motion, the boundary action becomes of the form 



S 



B 



K / d 3 x- 



u 2 aMu 



-K / dx + 6(x 



2tt 2tt ( 



-aM 



d 2 y 



d 2 k d 2 k' 
(2tt) 2 (2tt) 2 



,i(k-k')-y 



x J{ k 'a)(^i{uB,k' k p )J{k ll ) +F 2 (u B ,k'k IJ ,)J(k IJy ] 



(28) 



where J-"'s are given by the equations (14) and (15). 



Note that 6(x + ) represents our restriction on the region of time to be x + > 0, which 
is very natural for our setup from the metric point of view. The physical significances and 
effects of the time boundary and associated boundary condition are explained in detail in 

mi. 

To evaluate the momentum space correlation functions, we first evaluate the time integral 
for the range < x + < oo and define Q(w) = J dx + e lwx+ 8(x + ) ■ (ax + )~ aM . Then 



Q(w) = a 



-aM! 



W\ 



T(l - aM) -i cos( — - — )sign(w) + sm( 



(29) 



2 ' ° v ' v 2 
where w = u>' — u. These results are the same as those evaluated by Laplace transform. 
Thus the momentum space correlation function is given by 

^(uj\k')<p(oj,k)) = -2(27r)" 3 (5(A r/ - k) Q{J - u) F x {u B ,uj,k) , (30) 



where T\ is given by (17). We use the spatial translational invariance to get the delta 
function for k, but not for uj due to the lack of the time translation symmetry. 

Let us calculate the coordinate space correlation function as 



(cf)(xt,y 2 )<f)(xf,y 1 )) 

du' d 2 k' du d 2 k 



(31) 



2tt (2tt) 2 2tt (2tt) 2 
After carrying out the integral for the u/ first 

du' 



(a 2 x^x^") 2 (0(u/ ,k')(f)(uj,k)) . 



2tt 



(xt) a~ aM 



\—aM 



(32) 



the expression reduces to the Schrodinger case in the time independent part, whose calcula- 
tion is done similarly. Finally we obtain 



xt\ 2 f duj d k 



(0(4,y 2 )0(4,yO) = -20(4) (^J J — —e-^^e^-^T^u^) 

a Ad 

v L 3 M 1+U 0(4)004 - xt) fx+\ ^ ( M(y 2 -m) 2 \ 

■ ex P n , + +v ■ y 66 ) 



~~ I» ti2^u A z 2u (xt - xt) 2+u \xtJ V 2(4 - xj 

This result has exactly the same form as that of the aging case computed in |llj[T0 
We can similarly evaluate (00) correlation function as 

(0(4) 2/2)0(4 > yi)) 

aAd 

1 L 3 M 1+U 0(4)0(4-4) fx+\ 2 exp f M{y 2 -y^ 2 



T(y)TxV-^ B 2v (x+-x+) 2 +» V4/ V 2(4-4) 

f r(i/)' 



2(2+ -x; 



(34) 



This expression is our main result, which has all the essential properties related to the 
logarithmic extension of the aging holography. Let us comment about this coordinate space 
correlation function. 

• There are clear differences between the two functions (0(4) ^)0(4> $0) anc ^ (0(4' $2) 
0(4) The latter can be obtained simply by taking a derivative of the former with 
respect to the parameter v. This is clear from the structure of the differential operator. 
The same result can be obtained by a direct, yet lengthy, calculation outlined in the 
appendix £}A} 



The correlation function (0(4,2/2)0(4)2/1)) reveals an additional logarithmic depen- 
dence only of the time difference ln(4 — 4); besides some irrelevant constant factors. 
This is the main feature of the aging LCFT from our bulk construction. Interestingly, 
the spatial dependence y 2 — yi only shows up in the exponential part, which is respon- 
sible for the well defined decaying properties of the Schrodinger and aging correlation 
functions for both the large distance y 2 — y\^-oo and small time xt — 4 — regimes. 

We confirm the relation between the aging correlation functions and those of the 



4 We have a slightly different expression compared to [TT] using the identity of the T function. This form 
is slightly more convenient for the evaluation of Ti- 
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Schrodinger 



( 




) 



2 



(H x t>ih)4( x l>yi)) Aging 



(35) 



( 




) 



2 



(<t>{xti y2)4>( X l , Vl)) Aging 



(<f>(xt,y2)(f>{xf,yi))schr , 



(36) 



which is one of the important properties of aging observed in [10] [IT] . 

• We stress that the final form of the correlation function with logarithmic extension 
is fixed by the backgrounds motivated by string theory solutions. Schrodinger and 
aging backgrounds can be constructed directly from the AdS 5 x S 5 solution using the 
null Melvin twist [46j[47] for 7^0 and the light-cone reduction j3H]|13] for 7 = 0. 
Attempts to construct the corresponding dual field theories have been put forward, 
e.g. in 00] EE- 

4 Physics of two-time correlation and response functions 

We obtained the correlation functions and their logarithmic extensions of the Schrodinger 
case in £j2] and aging case in §[3} In this section, we describe the physical properties of 
them in detail and seek their applications to the KPZ universality class, motivated by its 
recent theoretical extension to the aging regime [Hj. Aging cases offer more general class 
of correlation functions with explicit breaking of time translational invariance compared to 
the Schrodinger case. Thus we concentrate on the aging phenomena and comment on the 
Schrodinger case briefly. 

Let us first comment about the physical identifications of (00) and (00). From the 
point of view of our detailed construction outlined in £|2j these are naturally identified as 
correlation functions of the logarithmic pair and 0. Yet, it is compelling to identify them 
as response functions by coupling the logarithmic pair (0, 0) with some external fields (h, h) 
as 8C = h ■ + h ■ 0. Starting from one point function (0) , the variations with respect to h 
and h will give the response functions (00) and (00), respectively. Similar identification in 
the field theory side is already adapted in the context of aging and its logarithmic extension 
[8 J [9J. Thus we would like to identify (00) and (00) as our response functions. In this section, 
we freely take the view them as correlation functions or response functions. We investigate 
the properties of (00) and (00) in this section and try to seek possible connections to KPZ 
class. 
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4.1 Correlation function 



The two-time correlation function of the geometric realizations with aging isometry is iden- 
tified and analyzed in the aging regime [TT], reproducing the characteristic features of its 
time dependence: slow power-law decay, broken time translation invariance and dynamical 
scaling [12] . Here we would like to analyze the full profile of the correlation function without 
the logarithmic extensions, motivated by the recent interest in the field theory side [8J[9J. 
As already noted above, the response function (00) has the same form as the correlation 
function obtained in [UJ. Thus we analyze its physical properties here together. 

Let us consider the two-time correlation function C(x2,xf) = (0(^2~, y\)<f>(xi, yi)) ob- 



tained by identifying y 2 = y\ from (33) 



aA4 
2 



C(x+, xt) = C (x+ - x+y 2 -" • 1^1 , (37) 
where Co = j^rx iL^t-iv ■ This depends on two different time scales, x\ and called 

\ / 7T^s XL g 

a waiting time, which is the time passed after the system is quenched or prepared. Now the 
system is disturbed with small perturbations at this waiting time xf, and subsequently 
measured at x\ which is called a response time. Explicit dependence on the waiting and 
response times in the two-time correlation function is a crucial difference of the aging system. 
The correlation functions of the Schrodinger invariant system or of the systems with time 
translational invariance, in general, only depends on x\ — xf, the time difference between the 
response time and the waiting time. These time dependent features manifest themselves in 
the term with aAi dependence. The correlation function without the logarithmic corrections 
is explicitly depicted in [11] for the aging regime, which is x\jx\ ^> 1 and x\,x\ ^> 1. 
There, it is demonstrated that the plots show the slow power law dynamics with explicit 
waiting time dependence, in addition to the excellent collapse in the aging regime revealing 
the universal features. 

Here we would like to explore the entire range of the waiting and response times motivated 
by recent interesting applications [8J [9 J . It is more transparent to adopt the parametrization, 
V = x t/ x ti un der which 

C(s, y) = C s-^y- 2 -^ (l - ^ ^ , (38) 



where s = x{ is the waiting time, 1 < y < oo and v = a/4 + m 2 L 2 — 7.M 2 . Note that it 
is only a function of the scaling variable y apart from the dependence of the waiting time 



s - 2 -r 
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4.2 Correlation function, roughness and exponents 



Growth phenomena have attracted much attention due to their wide applications in physics 
and beyond. Their examples include paper wetting, epidemics, flux lines in a supercon- 
ductor and fluid flow in porous media, see e.g. [4|[5j[lJ. We would like to point out two 
characteristic properties of these phenomena, which are relevant for the following discus- 
sions. First, this is intrinsically nonequilibrium process and breaks the detailed balance p], 
while the time independent treatments still provide much insight [!]. Simply the system is 
out of equilibrium due to the fact that noise or randomness is an essential part of its dy- 
namics. Second, the concept of the local scale invariance and associated universal exponents 
are useful descriptions, which are similar to the equilibrium critical phenomena. The second 
property is particularly important for our purpose. We consider the aging as well as the 
growth properties both for the equilibrium and for the out of equilibrium systems, using the 
Schrodinger and aging holography, respectively. 

Recently, Kardar-Parisi-Zhang type universality class in 1+1 dimensions is experimen- 
tally realized in the growing interfaces of the turbulent liquid crystals [6] j7] . The roughness 
of interfaces is quantified by their width w(l,t) defined as the standard deviation of the 
interface height h(x, t) over a length scale I at time t, or equivalently the radius R(x, t) 
along the circle due to the fact that the universality class is independent of the shapes of the 
interface [7]. The authors measure the interfacial width w(l,t) = (a/ {[h(x, t) — {h)i\ 2 )i) and 
the height-difference correlation function C{l,t) = ([h(x + l,t) — h(x,t)} 2 ) where (•••); and 
(■ • • ) denote the average over a segment of length I and all over the interface and ensem- 
bles, respectively [6]. Both w(l,t) and C(7,t) 1//2 are common quantities for characterizing 
the roughness]^] for which the so-called Family- Vicsek scaling [18] is expected to hold. The 
dynamical scaling property is 



with two characteristic exponents: the roughness exponent a and the growth exponent b. 
The dynamical exponent is given by z = a/b, and the cross over length scale is U ~ t x l z . 
For an infinite system, the correlation function behaves as C(l,t) 1 / 2 ~ t b at large times. 

In the holographic approach, the dynamical exponent z and the number of spatial direc- 
tions d are built-in. In our case, z = 2 and d = 2. Thus we can fix the roughness exponent a 
by determining the growth exponent b from the limiting behavior of the correlation function, 
C(l, t) 1 / 2 ~ t b , at large times for our infinite system. Once we take \y2— Vi\ — I < t = x^—x^, 

5 We concentrate on the correlation function due to the fact that both the w(l,t) and C^jt) 1 / 2 show 
similar scaling properties, and the latter is readily available for the holographic theories. 




for / /* 
for / > 



(39) 
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the correlation functions (19) and (33) behave as 



(4>(xt + t,l)(f>(xt,0)) Schr - 

{(f>(xt + t, l)4>(xl, 0)) Aging 



-2-v 



2-v-'- 



where we assume t 3> x± in the latter. Thus 

kschr = — 1 — , &Schr = 

b Aging = "1 ~ V jl - OlM / A , A Aging 



-2 



-2 - v - aM/2 



(40) 
(41) 



(42) 
(43) 



where the ratio between a and b is fixed, a/b = z = 2!r\ These relations show that the 
growth exponent bschr in the Schrodinger case is given by the parameter u, which measures 
the conformal weight of the dual scalar operator, while the dynamical exponent is fixed from 
the construction as z = 2. Similar story holds for the aging case, where the growth exponent 
b Aging is a function of ahA in the metric, as well as v. Note that there are interesting 
qualitative changes for a = or b = 0, which become clear below. 



v=2.1 



v=2.1 



C(y,s) 



itr-- 



C(y,s) 



Aging 



i» 4 



,2o y 




Figure 1: Left: Log-log plot of the correlation function C(y,s) vs. y given in (38 1, which shows a typical 
behavior of aging and is compared with the typical aging plot from numerical study, e.g. fig. 1 of [8]. Right: 
Log- log plot of C(y, s) vs. y — 1. 



With these discussions in mind, let us investigate our correlation functions in detail. First 
of all, it is natural to have v > 2 in our construction. For this range, we depict the log-log 
plot of the correlation functions for the Schrodinger and aging cases together in figure [T] 
We remind that the roughness w(l,t) shows the similar behaviors as the square root of the 
correlation function, C(/,t) 1 / 2 . In the left plot, we observe the typical aging behavior, which 

6 Surely, these exponents are not belong to the KPZ universality class (a = 1/2 and b = 1/3 in 1 + 1 
dimensions), which is realized in growing interfaces of the turbulent liquid crystals [6][7]- It is necessary 
to study the holographic theory with z = 3/2 to directly contact with the KPZ class. Our non-relativistic 
conformal case with z = 2 is more close to the Edwards- Wilkinson class 49J with the same dynamical 
exponent. 
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can be compared to fig. 1 of [8]. In the right plot, we put both the Schrodinger and aging 



correlation functions to see clear differences. Note that our aging correlation function (38) 



changes its behavior around ln[y — 1] ~ 0, while that of the Schrodinger does not. The latter 
represents the characteristic features of the system with translational invariance. 



B : v=-1.9999 



C(y,s) 



C(y,s) 

10"r 



A : v=-1.9 




r,y-i 




C(y,s) 




Figure 2: Log-log plot of the correlation function C(a 



We clearly observe the growing or decaying 



patterns. The top left panel (A) shows only decaying behaviors, valid for —2 < v 1 while the bottom right 
panel (D) shows only growing ones valid for v < -2 - sM. The bottom left panel (C) shows an interesting 
behavior, changing from growing to decaying, valid for —2 — 2^ < v < —2. Qualitative change is observed 
near v = — 2 as shown in panel B. 



It is interesting to extend our analysis to the range v < 2 from the phenomenological point 
of view. As mentioned above, there are qualitative changes at v — —2 and v = —2 — ^-. Let 
us concentrate on aging case, whose correlation function is fixed by our bulk geometry. It is 
interesting to see the qualitative different behaviors of the correlation functions by changing 
the values of parameters, v and <xM. The typical behaviors are depicted in figure [2j Their 
behavior for v > —2, given in panel A of figure |2j is similar to that of the v = 2.1 depicted 
in figure [TJ As v approaches the critical value —2, it becomes almost flat for small, positive 
ln[y — 1]. On the other hand, it shows the growing behavior for v < —2 



aM 
2 • 



14 



Surprisingly, the bottom left panel (C), which is valid for the range —2 



aM 



< V < 



shows a typical growing behavior for ln[y — 1] < 0, while decaying (aging) behavior for 
ln[y — 1] > 0. We can see these changes because we choose the horizontal axis as ln[y — 1]. 
If we choose the horizontal axis as ln[y] instead of ln[y — 1], we only see the aging behaviors 
because the entire region ln[y — 1] < shrinks to very small part around the origin of hx[y]. 



4.3 Response function 



Let us consider the two-time response function R^x^xf 



, which includes the loga- 



rithmic contribution obtained from (|34|) 

R{ x 2 ) ^1 ) 



■ L 2 x l ) 



.1 n 



2 



(44) 



^Y{v) ■ This is again the function of the waiting time and 



where R = 1/u + \n[Mu%/2] 
response time. See some relevant explanations on them near the equation (38). 

R(y,s) v=2.1 

4 x 10 63 r 



3 x 10 63 - 



2 x 10 63 - 



1 x 10 63 




0.000010 0.000015 0.000020 0.000030 



0.000050 0.000070 



0.00010$ 



-1 



Figure 3: Log- linear plot of the two-time response function with logarithmic extension (45) for v = 2.1 
with fixed aM. = 1. This parameter range of v naturally arises from our holographic model. 

To analyze R(x2 ,x^), it is convenient to introduce a scaling time variable y = x\jx\^ 
which is the ratio between the response and waiting times. Then 

R(s,y) = s- 2 -»y- 2 -"-^ (l - - Ins - Iny - ln[l - M , (45) 



where, once again, s = x± is the waiting time, 1 < y = x\ jx\ < oo and v = a/4 + m 2 L 2 — 
Typical behaviors of response function with logarithmic extension from our model are de- 
picted in figure [3j Similar behaviors of aging with logarithmic extension are noticed in 



15 



It is interesting to survey the general features of the response function R(s,y), for the 
different values of the parameter v as well as for the entire range of y, as a function of y. 
Before proceed further, note that our correlation function has explicit ln[l — \jy\ and lny 
dependence in addition to Ins as a consequence of the logarithmic extension. One might 
wonder the effect of lny which has the divergent behavior for large y. It turns out this 
depends crucially on the value of the parameter u, which is explained at the end of this 
section. 



B : v=-1.9999 



A : v=-1.9 




1 100 



C : v=-2.1 



R(y,s) 

120 000 r 




- 

io- 15 



R(y,s) 

4000 r 



D : v=-2.5 




i(r 5 o.oi io io 4 io' 



100 10* 



v y-i 



Figure 4: Plot of the two-time response function with logarithmic extension (45) for four different values 
of v with fixed aAA = I. The horizontal axis is y — 1. Each panel has three graphs with different waiting 
times. They are normalized to be unity at y = 10 10 to show the scaling behavior, and thus the effect of the 
factor s~ 2 ~ v is washed out. The remaining effects of the waiting time are from Ins in (45). The top left 



panel A is typical for v > — 2, which converges to for large y and diverges toward the negative infinity for 
small y. The panel B is for near v = —2, which shows the transition behavior of R(s,y). Panel C shows a 



typical behavior for —2 



< v < — 2, which nicely converges for both large and small y. Panel D is for 



aM 
2 



, which shows a diverging behavior for large y and thus is normalized at y = 10 instead 



There exist essentially four different types of qualitatively different behaviors, two from 
large y and two from small y. At large y, 1 — 1/?/ becomes a unity. Thus the behavior of 



16 



R(s,y) is controlled by the combination 



-2-u- 



2 my 



(46) 



which converges for u > —2 — and diverges otherwise. At small y — >■ 1, the behavior of 



R(s, y) is controlled by the combination 



l-l/!/)- 2 -ln[l-l/y] 



(47) 



which converges for v < —2 and diverges otherwise. From these behaviors, we conclude that 
the two-time correlation function converges both for small and large y for the range 



-2 - aM/2 < v < -2 . 
This case is depicted in the bottom left panel (C) of figure |4 



(48) 



Let us comment the effect of Iny, which would diverge for large y, as we mention above. 
Certainly, there exist a divergent behavior at large y for the range v < —2 — as depicted 



in the bottom right panel (D) of figure |4J For the other range, the effect of the lny does not 
produce qualitatively different behaviors in the response function. We investigate this issue 
by plotting the correlation function without the In y term 

-2-u 



Rl(s,v) = 8 



-2-Uy-2-u- 



C -lns-ln[l 



(49) 



against the two-time response function R(s,y) given in (45). The plot is depicted in figure 
[5| As expected, there are some visible differences for large y. But the differences do not 



A : v=-1.9 



B : v=-2.1 



R(y,s) 




100 10 4 



with Log y 
without Log y 



y-i 




with Log y 
without Log y 



Figure 5: Plots for the two-time response function (45) against Rl(s, y) without the lny term given in (49) 
for two different values of v. They are normalized by the value of (45) at y= 10 10 to show clear differences. 
The left panel is a typical case for v > —2, while the right panel is for —2 — <v< —2. For both cases, 
the effect of Iny is negligible for small y, while there are visible differences for large y. Even taking into 
account of all these differences, the qualitative features remain intact for v > — 2 — and thus physical 
applications would not be affected much. 
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change the qualitative features of the plot, and accordingly physical applications would not 
be affected much by the effect of \ny. 

With these observations, we go on to investigate the possibility that our system might 
be relevant for the physics of logarithmic extension of aging in the following section. 



4.4 Logarithmic extension of aging in KPZ class 



There have been various attempts to understand the logarithmic extension of the non- 
relativistic AdS/CFT [37J. Recently, the logarithmic corrections to the aging correlation 
and response functions are introduced to investigate their connections to the aging phe- 
nomena in field theory context, where their dynamical scaling is characterized by the aging 
exponents [8j[9j. Specifically the authors of [8 J study the aging phenomena of the surface 
growth processes described by (l+l)-dimensional Kardar-Parisi- Zhang equation [3]. The 
associated dynamical exponent is 3/2. 

While we are very interested in providing detailed analysis of our model, quantitative 
comparisons are not feasible due to several reasons at this point: our model is designed for 
2 + 1 dimensional dual field theory with dynamical exponent z = 2 exactly at the critical 
point Q We already presented our qualitative results with the logarithmic contributions in 
the previous section. Let us comment the similarities and differences compared with [8J[9J. 
Especially, the equation (4.3) of [9] (or similarly equation (10) of [8]) reads (using the same 
notation as in (45)) 

R(8,y) = 8- 1 ^y- x *t'(l-'^\ (50) 
x (ha - gi2,o£' ln[l - l -\ - l -M' 2 ln 2 [l - ^ - <? 2 i,o£' Hv ~ 1] + \m' 2 ^iv ~ 1]) , 

where the parameters #12,0, 921,0, fo come from the logarithmic extension in the field 
theory side. From the phenomenological input |9j "the parenthesis becomes essentially con- 
stant for sufficiently large y" the condition £' = is imposed to remove the last two terms. 



The other parameters in (50) are fixed to the data. Compared to our result (45), we can 



identify the various exponents 

1 aM \ R 
v = a — 1 , — — = 1 — a , (51) 

along with the condition a 1 = a, which is a restricted case. 

7 The case with z — 3/2 for arbitrary d + 1 dimensions has been considered in some detail [50]. There, 
correlation function without logarithmic extension is obtained for the special case, without full generality 
due to the technical difficulties. 
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There are also three notable differences in the parenthesis. First, we note the existence of 
the term Iny, which might spoil the large y behavior. As we discussed in the previous section 
the response function (45) converges at large y for the parameter range v > —2 " Vi 



U.3 



2 



Thus there exist large parameter spaces we can reliably describe the desired properties of 
the response function. Second, we also have an explicit dependence of the waiting time s. 
Even after the normalization with s 2+u C(s, y) to remove the dependence on the waiting time, 
there exist another source of the contribution on s, which is the important consequence of 
the broken time translational invariance and might provide interesting physical properties. 
Third, the correlation function (45) does not have ln 2 [y — 1] term. Seemingly this is a small 
correction and thus does not give qualitatively different features] 



We don't claim that the response function (45) quantitatively describes the scaling prop- 
erties of the growth phenomena for several reasons, including the different spatial dimensions, 
different dynamical exponents and different terms present. As far as we understand, there 
exist tensions between the experimentally realized growth and the theoretically investigated 
aging phenomena for the KPZ universality class. Our model provides several qualitatively 
distinctive behaviors of the correlation and response functions. In a particular case, the 
two-time correlation and response functions reveal an initial growth behavior as well as the 
slow power law decay, typical for aging, at later times for some parameter ranges. Further 
investigations along this line seem to be very interesting. In particular, studying the 1 + 1 
dimensional aging system with logarithmic extension for z = 3/2 would be very interesting 
for obvious reasons. 



5 Conclusion 

Physics out of equilibrium is challenging, yet very important. The aging system in holography 
as well as in field theory provides a simple framework describing time dependent, thus far 
from equilibrium, physical systems with local scaling symmetry. It has been known that 
corresponding correlation functions reveal a slow dynamics with power law decay, broken 
time translational invariance and dynamical scaling in the aging regime, where waiting time 
s, response time t and the scaling time y = t/s are large compared to the microscopic time 
scale of the system considering [2][TT]. 

1+1 dimensional Kardar-Parisi-Zhang (KPZ) universality class [3] was realized recently in 

8 In the context of the log CFT and its geometric realizations, an extension to the tri-critical gravity is 
shown to produce the In 2 [y — 1] contributions [5T] . Thus it is reasonable to expect that similar generalization 
would give the term in the case of Schrodinger holography, too. We thank to Y. S. Myung for the comment 
on this issue. 
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growing interfaces of turbulent liquid crystals, which is described by the dynamical exponent 
z = 3/2, growth exponent a = 1/2 and roughness exponent b = 1/3 [6] [TJ. Subsequently, 
aging properties of the same KPZ universality class in 1+1 dimension are studied in [S] with 
logarithmic extension of the correlation and response functions, where the critical exponents 
are further extended. These results serve as our motivation. We extensively study our 
holographic model to reveal the qualitative features of correlation and response functions 
and to seek the connection to the KPZ universality class. 

In this paper, we investigate the entire range of the scaling time y without assuming 
the aging regime, especially the early time profiles of the two-time correlation and response 
functions. We extensively study the typical patterns of them for the different parameter 
range with phenomenological point of view. In physical terms, the roughness of the grow- 
ing interface and the square root of the height-difference correlation function show similar 
behaviors in the KPZ universality class. There are three qualitatively distinct roughness 
behaviors in log-log plots summarized in figure [2j i) slowly decreasing in early times and 
continue decreasing faster in later times, ii) increasing in early times and change into de- 
creasing in later times and iii) fast increasing in early times and continue increasing slowly 
in later times. In particular, the second case shows an initial growing behavior, typical for 
the growth, followed by power law decaying behavior, typical for the aging for the parameter 
range —2 — 9 ^ < v < —2, which is unexpected and interesting. For the Schrodinger case, 
there are only two exclusive options, either growing or aging, in contrast. 

It will be interesting to construct the corresponding holographic model to check whether 
the holography can reproduce these critical exponents and provide more useful information 
on this universality class. To do so, it is required to extend our construction to non-conformal 
cases for z = 3/2inl + l dimensions, where analytic results are not available yet. 
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Appendix 



A Explicit evaluations of the logarithmic extension 

In this section, we show the details of the inverse Fourier transform of T%- The equivalent 
and simple derivations are presented in sections ^2] and §|3] by taking the derivative of T\ 
with respect to the parameter v. 



From (14), the non-trivial contribution of Ti can be expressed as 



* = (11 (<^(««») I, < 52 > 

where we expand the expression for small ub- The inverse Fourier transform (IFT) of q 2u in 
the first term is the same as T\. Thus we concentrate on IFT of q 2v In q 2 in the second term. 

By defining Q = q 2 = k 2 + 2iAiu, we compute IFT with respect to u as 



oo poo 

e iujt q 2u \nq 2 du = / e^Q v hxQduj 

-oo J — oo 

2 \ i>k 2 +i-oo 



2\M\ 



e A-^LZ e A^ Q ^ nQdQ - (54) 



For Ai > 0, a branch point is located at Qo = 0, while Uq = ^ is in the upper half plane. 
To get a nonvanishing contribution for positive t, we choose our contour to the upper half 
part in complex u plane. This is depicted in the left side in figure [6j Similarly, in the Q 
plane, the contour is left half plane of the line Re[Q] = k 2 and branch cut lies in the negative 
direction from the origin (Q = 0), which is depicted in the right part of the figure [6] 



From the contour integral in complex Q plane, the nonvanishing contributions come from 
two terms 

° + +. exp (2^) Q u ^Q d Q + f Q °° " ex P (2^) Q u ^Q d Q ■ (55) 

In the first term we identify Q = \Q\e l7T , while Q = \Q\e~ l7T in the second term. Thus we get 
2ni cos(7rz/) J exp (^~~^J^Q^j Q u dQ + 2ism(nu) J exp f—^-^Q^j Q" In Q dQ , (56) 
which can be evaluated to be the following using the definition of the T function 



r> + r 



2MY +1 L . , , ~ 2M o . . . .r(i/+l) 



2ni cos(7ri/) + 2i sin(iih') In h 2i sininu 



t v Tv+1 



(57) 
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Note that it is independent of the momentum k. Thus the integration for k can be done 



with the coefficient of equation (54) 



<^+k 2X2 ) e -^{kl+kl) dkidk2 = V^L^-fM^D . (58) 
As a result, the final form of the inverse Fourier transformation is simply given by 

2 2 

/ 7 — w> — e lk v +lu}X q 2u lnq 2 = — T(u + l)sm(uir) —— . (59) 

J (2tt) 2 2tt h h g7i2M dv i v J k J y x+ j j v j 

Thus we show that the correlation function (0(x^, jj2)(j)(xf, yi)) can be directly computed 
using the inverse Fourier transformation to confirm the same result obtained by the v deriva- 
tive 

M*+ ft)> = | ^^e*&^W*^r a {u B ,u,,%) (60) 

#(Ax+)L 3 _ma£ ( nu (2M \ u+2 /u B \ 2v \ 

exp 2Ax+ (61) 



8n 2 uu%M du \ T(u) \Ax 

~{4>^t,m^i,yi))- (62) 
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